Spin glasses in a field: Three and four dimensions as seen from one space dimension 
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We study the existence of a line of transitions of an Ising spin glass in a magnetic field — known 
as the de Almeida- Thouless line — using one-dimensional power-law diluted Ising spin-glass models. 
We choose the power-law exponent to have values that approximately correspond to three- and four- 
dimensional nearest-neighbor systems and perform a detailed finite-size scaling analysis of the data 
for large linear system sizes, using both a new approach proposed recently [Phys. Rev. Lett. 103, 
267201 (2009)], as well as traditional approaches. Our results for the model corresponding to a 
three-dimensional system are consistent with there being no de Almeida-Thouless line, although the 
new finite-size scaling approach does not rule one out. For the model corresponding to four space 
dimensions, the new and traditional finite-size scaling analyses give conflicting results, indicating 
the need for a better understanding of finite-size scaling of spin glasses in a magnetic field. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

One of the most striking predictions of the mean- 
field theory of Ising spin glasses, taken to be the exact 
solution 1 - of the infinite-range Sherrington-Kirkpatrick 
(SK)2 model, is the existence of a line of transitions, 
known as the de Almeida-Thouless 3 - (AT) line in the pres- 
ence of a magnetic field. This line of transitions separates 
a high-temperature region where the description of the 
model is quite simple, just involving a single order pa- 
rameter, from a low-temperature region where there is 
"replica symmetry breaking" (RSB) in which the system 
has an infinite number of order parameters characterized 
by a function^ 

The question of whether RSB applies to realistic short- 
range spin glasses remains controversial. According to 
the RSB picture, real spin glasses behave rather simi- 
larly to the SK model and so have an AT line. However, 
according to the phenomenological "droplet" pictures- 
there is no AT line and the zero-field transition is rounded 
out by any nonzero magnetic field in finite-dimensional 
short-range systems. 

It is convenient for simulations that a static property, 
the spin-glass susceptibility x SG , diverges at the transi- 
tion. This quantity is the inverse of the eigenvalue of the 
stability matrix* 3 - found by de Almeida and Thouless and 
is given by the zero- wave- vector k — limit of Xsg s ) 
defined in Eq. (O below. Because x SG can be computed 
in simulations directly^ one might imagine that it would 
be straightforward to decide if the AT line occurs in, say, 
a three-dimensional spin glass. However, there is still no 
consensus on this issue because there seem to be quite 
large corrections to finite-size scaling (FSS). The purpose 
of this paper is to investigate different methods that have 
been proposed to perform FSS to see if there is an AT 



line in three and in four space dimensions. 

In fact, rather than to study short-range models^— 
we find it convenient to study a one-dimensional model 
with long-range interactions which is taken as a proxy for 
a short-range model.— ~— The interactions of the long- 
range model fall off with a power a of the distance and 
varying the power effectively corresponds to varying the 
space dimension of the corresponding short-range model. 
Here we study long-range models which are proxies for 
short-range models in three and four space dimensions. 

The outline of this paper is as follows. Section [TTJ de- 
scribes the model, the quantities we calculate and details 
of the numerical simulations. Section Mil describes the 
results, and Sec. IIVI gives our conclusions. 



II. MODEL, OBSERVABLES & NUMERICAL 
DETAILS 



A. Model 

We study a variation of the model introduced in 
Ref. Hi], which is given by the Hamiltonian 



£ ij J ij $i 



- 2J hi Si 



(1) 



In Eq. ([T]), Si = ±1 arc Ising spins placed on a ring 
of length L to enforce periodic boundary conditions in 
a natural way* 1 ^ The interactions are chosen from a 
Gaussian distribution with zero mean and standard de- 

i 3 takes values and 
t7a a of taking the value 



viation unity. The dilution matrix e 
1 , and has the probability 



1, where r,j = (L/tt) sin(7r|? — j\/L) is the geometric dis- 
tance between the spins. To prevent the probability of 



2 



placing a bond between two spins being larger than 1, a 
short-distance cutoff is applied and, thus, we take 



Pij = 1 - exp{-A/rl?) 



(2) 



The constant A is determined by the requirement that 
the mean coordination number, z av , takes a specified 
value 



L-l 

Yl Pi 

i=l 



(3) 



and we set z av = 6. The values of A are given in Table 
UJ The site-dependent random fields hi are chosen from a 
Gaussian distribution with zero mean and standard de- 
viation H . 

By tuning the exponent a in Eq. @ one can change 
the universality class of the model in Eq. (Q]) from the 
infinite-range to the short-range universality case. For 
< a < 1/2 the model is in the infinite-range universal- 
ity clas o 24 i 25 and, in particular, a = corresponds to the 
Viana-Bray models For 1/2 < a < 2/3 the model de- 
scribes a mean-field, long-range spin glass,— correspond- 
ing to a short-range model with a space dimension above 
the upper critical dimension, i.e., d > d u = 6 i 27 ' 28 For 
2/3 < a < 1 the model has non- mean- field critical be- 
havior with a finite transition temperature T c , while for 
a > 1, the transition temperature is zeroi^ 7 - The value of 
a in the long-range one-dimensional model corresponds 
roughly to an effective space dimension d in a short-range 
model via the relation 1 ^— 



d = 



2 ~ V(d) 
2a - 1 



(4) 



where 77(d) is the critical exponent 77 for the short-range 
model, which is zero in the mean-field regime. Here 
we are interested in values of a that correspond to 
three- (3D) and four-dimensional (4D) short-range sys- 
tems. Because rj(d = 3) = — 0.384(9) p£ three dimen- 
sions corresponds to ad=3 — 0.896, and since r)(d = 
4) = — 0.275(25)^ four space dimensions corresponds 
to a d=i ~ 0.784. 



B. Observables 

To determine whether a spin-glass state exists in a 
magnetic field, we study the wave- vector-dependent spin- 
glass susceptibility defined by 



(s,s s ) T - [SMS,), 



ik (i-j) 



(5) 

where (• • • )t denotes a thermal average and [• • -] av an 
average over the disorder. Each thermal average is ob- 
tained from a separate spin replica, i.e., we simulate four 
copies with the same disorder but different Markov chains 
at each temperature. As discussed in Sec. U x S g (ty i s an 



appropriate quantity to study because x S g = Xsg = 0) 
diverges on the AT line. 

It is also convenient to extract from the spin-glass sus- 
ceptibility a correlation length, which is usually defined 

^ ,.14, 18, 31, 32 



2sin(fci/2) 



Xso(O) 

XsgOO 



l/(2a-l) 



(0) 



where k\ = 2n/L is the smallest nonzero wave vector. 
Because we work in the non-mean-field regime, standard 
finite-size scaling (FSS) applies, i.e., 



£l/L = X\I}I v {T -T c ) 



(7) 



The importance of £l/L is that it is dimcnsionless. As 
such, data for £,l/L for different system sizes L cross at 
the transition temperature T = T c if corrections to FSS 
are unimportant [see, for example, Fig. HJa)]. This is a 
particularly convenient way to locate T c . 

For the one-dimensional model the critical exponent 
?7lr satisfies the exact relation 2 — ?7lr = 2a — 1, so it is 
also convenient to study the finite-size scaling of a second 
scale-invariant quantity, namely 



AV 



JL 2 -^ = XsJL 2 "- 1 = C[L^(T - T c )] 



(8) 
for 



Again, if FSS applies without corrections, data 
X SG /L 2cr_1 for different system sizes cross at T c . 

Both the correlation length in Eq. (O and the spin- 
glass susceptibility in Eq. (EI involve k = fluctua- 
tions. Recently, Refs. [13, [23l and l34l have argued that 
one should avoid data at k = for spin glasses in the 
presence of a magnetic field on the grounds that there 
are large corrections to FSS. We discuss this in detail be- 
low and for now just present the new proposed quantities 
to be measure d 22 ' 23 ' 34 that avoid k = fluctuations. We 
shall use the term "modified" FSS analysis to denote the 
use of these quantities. Below we compare the results of 
this modified analysis with results obtained from Eqs. ([7]) 
and (E]), which we denote the "standard" FSS approach. 

At long wavelength one expects 



Xsaik)- 1 ^XsaiOr'+By + Cy 2 



(9) 



with y = fc 2cr_1 , which is the generalization of the 
Ornstein-Zernicke equation to long-range interactions, so 
one can calculate X SG (0) 1 indirectly by fitting data for 
nonzero k to this form. For L — > 00, this extrapolated 
value, Xsg(^ — ^ 0) i vanishes at and below the transi- 
tion temperature. Interestingly, one find s 22 ' 23 that the 
extrapolated value goes through zero even for a finite 
system, at a temperature T*(L) which tends to T c for 
L — > 00. This means T c = lim£_>oo T*(L). We find quite 
strong corrections to the asymptotic result (as found also 
for the Sherrington-Kirkpatrick model in Ref. l35l ) , so we 
include the leading correction to scaling by fitting the 
data to 



T*(L) = T c 



A 



(10) 
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FIG. 1: (Color online) Data used in a standard FSS for a — 0.784, which is a proxy for four dimensions. The left column is 
for H = and the right column for H = 0.1. The top row shows the finite-size correlation length divided by L as a function 
of temperature T for different system sizes. The bottom row shows results for the scale-invariant spin-glass susceptibility, 
X SG /i 2<T_1 - The intersections seen in the left-hand column indicate a transition in zero field. The intersection temperatures 
(shown by the dashed vertical line) are T c ~ 1.34 from £l/L data in (a), and T c ~ 1.36 from xsg / 'L 2a ~ 1 data in (c). The small 
difference is likely due to sensitive corrections to scaling in the susceptibility (see Ref. [3^ for details). By contrast, the lack of 
intersections in the right-hand column indicates no transition in a magnetic field. 



where A is a correction to scaling exponent and A is the 
amplitude of the correction. 

This method requires several nonzero wave vectors and 
so is particularly suitable for one-dimensional models as 
studied here, because one can simulate very large linear 
sizes for these. For a short-range model in 4D, where 
the number of wave vectors is more limited, the authors 
of Ref. [34] propose another quantity. To motivate this 
quantity, we note as stated above, that £l / L is particu- 
larly useful because it is dimensionless. From Eq. ([6]) we 
see that the crucial quantity is the ratio Xsg(0)/Xsg(^i); 
which is also dimensionless. Massaging this quantity to 
obtain another dimensionless quantity ^ /L according to 
Eq. © is actually not essential. Therefore, a related 
quantity which does not involve k = can be defined via 



i?12 



Xs G ( fc l) 
Xso( fc 2) 



(11) 



where fc 2 = Att/L is the second smallest nonzero wave 
vector. Because it is dimensionless, R12 has the same 



FSS form as £l / L shown in Eq. ([7]) . Consequently, curves 
of i?i2 for different system sizes L should intersect at T c 
if corrections to FSS are unimportant. 



C. Numerical Details 

The simulations are done using the parallel tempering 
(exchange) Monte Carlo method^^ Simulation param- 
eters are listed in Table HI Equilibration is tested using 
the method developed in Ref. [ID [Eq. (8)]: The energy 
per spin is computed directly, as well function of a 

spin correlator. Both have to agree if a system is in ther- 
mal equilibrium. Starting from a random configuration, 
the directly computed energy is typically overestimated, 
while the energy computed from the correlator is under- 
estimated. Only when both agree (on average) is the 
system in thermal equilibrium. We thus perform a loga- 
rithmic binning of the data and requiring that both the 
energy per spin and the energy computer from the corre- 



4 



0.02 



0.01 - 



0.00 



-0.01 



1 1 1 


l | l l l l | l 


i i | i i i 


1 1 1 1 1 1 1 1 V 1 - 




h-ah L = 256 






; 


i-n L= 512 




if ''■ 




L = 1024 








* L = 2048 








h»h L = 4096 








o = 0.784 








H = 0.000 






: (a) i^-i^VV 

" ■ i ■ ■ j-<f ■ i 1 i - i i i i i i/ 





1.00 1.10 1.20 1.30 1.40 1.50 1.60 
T 



1.70 



1.40 



1.30 



_ 1 1 1 


i | i i i i | i i i i | i i i 


— 1 — | — 1 — 1 — 1 — 1 — | — 1 — 1 — 1 — 

a = 0.784 \ 








o 
o 
o 

ii 




^ L= 256 








L= 512 








h^. L = 1 024 








«h L = 2048 








L = 4096 






! (c) 

ijii 


— 1 1 1 1 1 1 1 1 1 1 1 1 1 1 L 





1.00 



1.10 



1.20 



1.30 
T 



1.40 



1.50 



1.60 



0.10 
0.05 

T 

a 

* w o.oo 

-0.05 



, , , I , , 


1,111,111 


i ■ ■ ■ i ■ ■ 'jn 


\ a = 0.784 




>-»h L= 128 




- H = 0.100 




L = 256 








* L= 512 


Jr/ ' 






>-»h L = 1024 








* L = 2048 




- 




h»h L = 4096 
















- (b) 












W i 1 i 1 i i 


S\ 1 i i M 1 i i i 





1.70 



1.30 
1.20 



1.10 

0. 



60 



0.80 



1.00 



1.20 
T 



1.40 



1.60 



0.60 0.80 1.00 1.20 1.40 1.60 1.80 
T 





-i — i — i — , — i — i — 


1 — , — 1 — 1 — 1 — | — 1 — 1 — 1 — , — 1 — 1 — 1 — £ 






a = 0.784 j 






H = 0.100 i 










L= 128 






* L= 256 






h»h L= 512 






i-«h L = 1 024 






* L = 2048 




1 (d) 


i-*-" L = 4096 








: , . , 1 


■ i i I i i 





1.80 



FIG. 2: (Color online) Data used in the modified FSS analysis for a — 0.784. As in Fig.[TJ the left column is for H = and the 
right column for H — 0.1. The top row shows the inverse spin-glass susceptibility XsgC — ^ 0) _1 extrapolated from results for 
k 7^ 0. The bottom row shows data for the dimensionless ratio R12 defined in Eq. (|lip . The dashed line for Xsg(^ — * 0) _1 [top 
row, panels (a) and (b)] indicates the extrapolated transition temperature T c determined from the fits in Fig. [3] The R12 data 
for H — [panel (c)] show intersections but the trend of intersection temperatures is not smooth for this range of sizes. The 
dotted vertical line shows the transition temperature estimated from the intersections of Xsg/^ 2<T_1 m Fig- Q] (Tc = 1.36) and 
so is just a guide to the eye. The R12 data for H = 0.1 [panel (d)] show no sign of a transition in the range of temperatures 
studied. 



lator agree for at least the last two logarithmic bins. The 
reason we wait for two additional logarithmic bins is be- 
cause the equality only holds on average. By being more 
conservative with the equilibration times we ensure that 
the bulk of the samples are in thermal equilibrium^ In 
addition, we verify that all other observables are inde- 
pendent of Monte Carlo time for at least these last two 
bins. 



III. RESULTS 

A. a = 0.784 (four space dimensions) 

When a = 0.784 the one-dimensional model is a proxy 
for a short-range model in four space dimensions. Data 
f° r £l/L an d Xsg/^ 2<t ^ u sed within the standard FSS, 
are shown in Fig. Q] for H = (left column) and H = 0.1 
(right column). The zero- field data for the scaled x SG in 
Fig. [ljc) show clear intersections indicating a transition 



at T c ~ 1.36, while the data for in Fig. Ufa) show 

clear intersections at T c ~ 1.34. This (small) difference is 
presumably due to corrections to scaling. Because there 
is no doubt that there is a zero-field transition for the 
models studied in this paper (see e.g., Ref. I20I ) we have 
not carried out a precise estimate of the value of T c in 
zero field. 

In contrast to the zero-field data, the data in a field 
of H = 0.1, Figs. QJb) andQJd), show no intersections, 
indicating the absence of a transition, at least for this 
range of temperatures and field. 

Data for the extrapolated value of x SG (k — > 0) _1 fitted 
according to Eq. © and R12 defined via Eq. (fTTj) are 
shown in Fig. [5] and used in the modified FSS. The zero- 
field data for R12 in Fig. [5{c) show intersections but the 
intersection temperatures do not vary monotonically for 
this range of sizes. The dotted vertical line in Fig. [2jc) 
corresponds to T c from the scaled x SG data in Fig. (O(c), 
and so is only a guide to the eye. The data for R12 in a 
field (H = 0.1) in Fig. [^d) have no intersections, i.e., no 
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FIG. 3: (Color online) Values of T*(L), the temperature where Xsg(^ — ^ 0) 
H = (left) and H = 0.1 (right) for a = 0.784. The data are fitted to Eq. 

In both cases a finite value of T c (the intercept) is found. For H = the goodness-of-fit paramete r 39 ' 40 Q is not very small 
indicating a satisfactory fit. For H = 0.1 the value of Q is smaller because the point for the largest size is well below the fit. 



1 goes through zero, for different system sizes for 
|10|) with T c , A, and A as adjustable parameters. 



transition is visible for the range of temperatures studied. 

For each system size, the temperature where the data 
for x S c(k ~* 0) — 1 snown m Figs. HJa) and HJb) goes 
through zero is referred to as T*(L). This value is ob- 
tained by fitting the data to a cubic polynomial (using the 
seven data points nearest to zero) , and error bars are ob- 
tained from a bootstrap analysis. The results are shown 
in Fig. [3] The thermodynamic transition temperature T c 
is then given by limL-^oo T*(L), which we estimate by 
fitting to Eq. (jTU)) . The fits, obtained by adjusting T c , 
A, and A, are shown in the figures (dashed lines). The 
data indicate a finite value of T c both in a field and in 
zero field. The central values for T c are shown by the 
dashed vertical lines in Figs.^ta) and[U[b). The value of 
X 2 /-/Vdof, where -/Vdof is the number of degrees of free- 
dom and the goodness-of-fit paramete r 39 ' 40 Q, indicate a 
good fit in the case of H = (Q = 0.50). For H = 0.1, 
the value of Q is smaller (Q = 0.014) because the point 
for the largest size is well below the fit. 

Our analysis to ascertain whether there is a finite T c 
includes four sets of data: £l/L, Xsg/ l L 2a ^ 1 , Xsc(k ~~ ^ 
0)-\ and R 12 . In zero field they all clearly show a fi- 
nite transition temperature, in agreement with the work 
of Banos et al^, who studied almost the same model. 
However, in a field of H = 0.1 there is an inconsistency: 
Three of the four measures (£,l/L, x sg /L 2ct_1 , and R12) 
show no sign of a transition. In contrast, the value of T c 
from results for x SG (fc — > 0) _1 does appear to be nonzero 
(T c = 0.91 ±0.07), see Fig.^b). We note, however, that 
the error bar is large and the last data point being below 
the fit may possibly indicate a downward trend at larger 
sizes. 

This discrepancy highlights the need to better under- 
stand FSS in spin glasses in a magnetic field. We shall 
come back to this question in Sec. IIV1 However, already 
we note that because the x SG — ^ O)^ 1 data are the only 
indicator for a finite T c in a field, we feel we should view 
results for this quantity with caution. 



B. a = 0.896 (three space dimensions) 

The one-dimensional long-range model with a = 0.896 
is a proxy for a short-range spin glass in three space di- 
mensions. The data used in the standard FSS analysis 
are shown in Fig. 01 Again, the left column is for H = 
and the right column is for H = 0.1. The zero- field re- 
sults for the scaled x SG m Fig- H{c) show a clear transi- 
tion. The zero-field data for £l/L in Fig- Hla) are less 
clear cut because there is little splaying of the data at low 
temperatures. However, the temperature where the data 
merge for two neighboring sizes increases as the system 
size increases. Similar results were obtained by Banos et 
al£2- for almost the same model, although they were able 
to study sizes of up to L = 8192 which do show a clear 
intersection with the data for L = 4096 (see Fig. 15 in 
their paper). Performing a detailed FSS analysis, Banos 
et al£Q showed that all their data are consistent with a 
finite value of T c . Because our data for do not in 

itself convincingly locate the transition temperature, the 
dotted line in Fig. B|a) shows (as a guide to the eye) the 
location of T c as determined from the scaled x S g data in 
Fig. He). 

In a field, the data for £l/L in Fig.@Jb) and the scaled 
spin-glass susceptibility in Fig.[4|d) show no intersections 
and, hence, indicate that there is no transition for this 
range of temperature and field. 

Data for the modified FSS analysis are shown in Fig. [5] 
The results for R12 are shown in the bottom row. In zero 
field [see Fig. [Sl[c)] there are clear intersections, although 
these do not vary monotonically for the range of sizes 
studied. However, the data are consistent with the value 
T c ~ 0.795 obtained from the scaled x SG data in Fig.fJJc), 
and this value is indicated as a guide to the eye by the 
dotted vertical line in Fig. EJc). The data for R12 in a 
field H — 0.1 [Fig. Hfd)] do not show clear evidence for a 
transition in the range of temperatures studied. 

The top row of Fig. [S] shows results for Xsg(^ — ^ O)^ 1 . 
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FIG. 4: (Color online) Data used in the standard FSS for a = 0.896, which is a proxy for three space dimensions. The left 
column is for H — and the right column for H = 0.1. The top row shows the finite-size correlation length divided by L 
as a function of T for different system sizes. The bottom row shows results for the scale-invariant spin-glass susceptibility, 
X SQ I . The intersections in the data for Xsg / l L 2a ~ 1 for H — [panel (c)] indicate a transition at T c ~ 0.795, which is 
shown by a dashed vertical line. The zero-field results for £l/L [panel (a)] show stronger finite-size effects and the evidence for 
intersections is not so clear. The same result was found recently by Banos et a/.—, who were able to study somewhat larger 
sizes for almost the same model. Their FSS analysis of all their data supports a zero-field transition. In panel (a) the dotted 
vertical line indicates the transition temperature found for the data in panel (c) and so is just a guide to the eye. The data 
in panel (a) are consistent with intersections tending to this value for L — > oo. By contrast, the lack of intersections in the 
right-hand column indicates no transition in a magnetic field. 



The temperatures where X SG (fc — > 0) 1 = are plotted 
and fitted according to Eq. (fTU)) in Fig. [6j The values of 
Q, 0.37 for H = and 0.91 for H = 0.1, indicate a good 
fit. For H = the result is T c = 0.81 ± 0.03 and the 
central value is indicated by the dashed vertical line in 
Fig. [3a). This value is consistent with the value 0.795 
from the data for the scaled x SG in Fig. |3Jc) ■ For H = 
0.1, the values of T*(L) shown in Fig. &[b) have a very 
strong size dependence. The figure also shows the values 
of the parameters obtained by fitting the data to Eq. (TTU1) . 
In particular, we find T c = —0.37 ± 1.15, indicating that 
the optimal T c is negative but the error bar is very large. 
This large error bar requires more discussion, which we 
now give. 

First, we note that a log- log plot of the data in Fig. [7] 
indicates that the data are compatible with T c — 0. In 



addition, we perform the following analysis: For each 
system size we construct 200 bootstrap data sets and 
estimate T*(L) for each of them. There is a huge scatter 
in the estimates from the largest size L = 1536. Thus 
we ignore this size in the analysis. For L — 784, 9 of 
the 200 bootstrap data sets do not yield a temperature 
where xsc(fc — > 0) _1 vanishes. Hence, we consider 191 
bootstrap data sets and fit each of them according to 
Eq. HQ]). 

Figure [5] shows the resulting cumulative distribution 
of transition temperatures, i.e., the probability that the 
transition temperature is less than the stated value. We 
find that 34 of the 191 data sets do not have a minimum 
in \ 2 ': rather, considering x 2 as a function of T c while 
optimizing with respect to the other fit parameters, x 2 
decreases monotonically as T c — > — oo while A — > in this 
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FIG. 5: (Color online) Data used in the modified FSS analysis for a = 0.896. As in Fig. [4] the left column is for H — and 
the right column for H — 0.1. The top row shows the inverse spin-glass susceptibility Xsg(^ — * 0) extrapolated from results for 
k 0. The bottom row shows data for the dimensionless ratio R12 defined in Eq. (|ll|l . The dashed vertical line in panel (a) 
marks the value for T c obtained by extrapolating the data for T*(L) [where Xsg(^ ~* vanishes] according to Eq. (|10|l ; see 
Fig. E] The result is T c — 0.81(3), in agreement with the value of 0.795 from the intersections of the scaled Xsg m Fig- HJc). 
The data for R12 in zero field in panel (c) show intersections, although the intersection temperatures do not vary smoothly. 
The dotted vertical line shows T c obtained from the scaled x S g ano - ^ s J us ^ a guide to the eye. The results for R12 for H = 0.1 
in panel (d) do not provide clear evidence for a transition. 
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limit. The estimate of T c from the global fit (T c = —0.37, indicated by the dashed vertical line) agrees well with the 
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TABLE I: Parameters of the simulations for different field 
strengths H and exponents a. N sa is the number of sam- 
ples, iVsw is the total number of Monte Carlo sweeps, T m i n is 
the lowest temperature simulated, and Nt is the number of 
temperatures used in the parallel tempering method for each 
system size L. The last column shows the parameter A in 
Eq. @ for z av = 6 neighbors. 
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FIG. 7: (Color online) Log-log plot of the data for T*(L) for 
a — 0.896 and H — 0.1 assuming T c — 0. The fit works very 
well according to the goodness-of-fit parameter Q, which is 
0.81. 
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median (50th percentile) of the bootstrap estimates. The 
median is indicated by a horizontal dashed line. Also in- 
dicated by horizontal dashed lines are the 16th and 84th 
percentiles, which would correspond to one standard de- 
viation if the distribution of T c 's were Gaussian, which 
is clearly not the case here. Only 30% of the bootstrap 
fits have a positive T c . We therefore conclude that a pos- 
itive T c is somewhat unlikely but cannot be completely 
excluded by the data for T*{L). 

To conclude this subsection, all data are consistent 
with there being no AT line at H = 0.1 for a = 0.896. 
The results for T*(L) obtained from the vanishing of 
X SG (fc — > 0) _1 do not exclude a finite T c , but this possi- 
bility docs not seem to be supported by the rest of the 
data. 



IV. CONCLUSIONS 

We have presented results of simulations of one- 
dimensional spin-glass models in a magnetic field which 
are proxies for short-range models in three and four space 



FIG. 8: (Color online) Cumulative distribution of the tran- 
sition temperature T c obtained from a bootstrap analysis of 
the data for a — 0.896 and H = 0.1. The estimate of T c 
from the global fit (T c = —0.37, indicated by the dashed ver- 
tical line) agrees well with the median (50th percentile) of the 
bootstrap estimates, which is indicated by a horizontal dashed 
line. Also indicated by horizontal dashed lines are the 16th 
and 84th percentiles which would correspond to one standard 
deviation if the distribution of T c 's were Gaussian, which is 
clearly not the case. Only 30% of the bootstrap fits have a 
positive T c . 

dimensions. We have analyzed the results using both a 
traditional FSS approach (which uses k — data) and 
a recently proposed modified FSS approach which uses 
only k > data. 

For the model which is a proxy for a 3D system, all our 
results are consistent with there being no transition in a 
magnetic field, at least for the range of fields and temper- 
atures that we can study. The results for T*(L), obtained 
from the modified FSS analysis, are just compatible with 
a finite T c but the other data are only compatible with 
the absence of a transition, at least assuming that the 
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data are in the asymptotic scaling region. 

For the model which is a proxy for a 4D system, three 
of the four sets of data indicate the absence of a transition 
in a field {^l/L, x S g/^ 2<t_1 ' anc ^ ^12)1 while that for 
T*(L) gives a satisfactory fit indicating T c > 0. This 
contradiction indicates that at least some of the data 
cannot be in the asymptotic scaling region. 

It is, therefore, crucial to understand whether it is bet- 
ter to use k = data in the analysis as in the stan- 
dard approach or to exclude that data as in the modified 
approach. 22 ' 23 ^ References I22I [23l. and [H argue that the 
k = data have strong corrections to FSS. On the other 
hand, the divergence occurs at k = and normally one 
uses divergent quantities in FSS because these should 
show the asymptotic FSS behavior for the smaller sys- 
tem sizes. It should also be noted that Ref. [ID found 
good agreement for the location of the AT line for the 
spin glass on a random graph, i.e., the Viana-Bray model 
(which corresponds to the a = limit of the present 
model), using the scaled x S gj i- e -i the k = fluctuations. 

Finally, an alternate interpretation of the results can 
be done using droplet scaling arguments^— The size £d 
of the droplets within this picture can be estimated by 
equating the domain- wall energy required to create them, 
~ J£q, to the energy which can be gained from flipping 

a droplet of size £d in the field, ~ H^ 2 , where d is 
the space dimension. For the long-range models studied 
here, 9 = 1 — <jH Thus, the droplet size is of order 3320 
for a = 0.784 and is of order 335 at a = 0.896, when 
H/J = 0.1, the ratio used in this study. We have only one 
data point, that for L = 4096 in Fig. [5Jb), greater than 
3320. Interestingly, it is the data points at the smaller 
values of L which point to a finite value of T c . The last 



point at system size 4096 lies well below the fitted curve. 
This suggests that had we been able to obtain data for 
a range of system sizes significantly greater than 4096, it 
might have been possible to obtain results in the analog of 
four dimensions like that displayed in Fig.[7]for the analog 
of three dimensions, where the last five data points are 
all greater than the estimated droplet size there of 335 
and the extrapolated value of T c is zero. 

Ideally, one would determine which set of data is in 
the asymptotic scaling regime by simulating larger sys- 
tem sizes. However, because the present study involved 
a rather substantial amount of CPU time, this is not fea- 
sible for us at present. Based on the data shown, the 
balance of the evidence is that there is no AT line in 
the one-dimensional models which are proxies for three 
and four dimensional short-range spin glasses. However, 
a deeper insight into corrections to FSS in spin glasses is 
needed to confirm this conclusion. 
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